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Abstract Almost 50 years after radio pulsars were discovered in 1967, our understanding 
of these objects remains incomplete. On the one hand, within a few years it became clear 
that neutron star rotation gives rise to the extremely stable sequence of radio pulses, that the 
kinetic energy of rotation provides the reservoir of energy, and that electromagnetic fields 
are the braking mechanism. On the other hand, no consensus regarding the mechanism of 
coherent radio emission or the conversion of electromagnetic energy to particle energy yet 
exists. In this review, we report on three aspects of pulsar structure that have seen recent 
progress: the self-consistent theory of the magnetosphere of an oblique magnetic rotator; 
the location, geometry, and optics of radio emission; and evolution of the angle between 
spin and magnetic axes. These allow us to take the next step in understanding the physical 
nature of the pulsar activity. 
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1 Introduction 


Radio pulsars are the archetypal observed neutron stars. Their discovery at the end of the 
1960s ( |Hewish et al.|[T968| l was definitely one of the major astrophysical events of the 
20th century. Their discovery confirmed the theoretical prediction of neutron stars in the 
1930s l |Landau|ri932[ [Baade & Zwicky||1934| l. Neutron stars have mass M of about 1.2- 
2.0 Mq, near the Chandrasekhar mass limit 1.4 Mq; and radius R of only 10-15 km. They 
result from the collapse of typical massive stars in the final stage of their evolution ( |Shapiro| 
|& Teukolsky ||1985|l; or from white dwarfs, when accretion from a companion star pushes 
them over the Chandrasekhar limit ([Whelan & Iben|1973[[Bailyn & Grindlay|1990[|Nomoto| 


|& Kondojl 199 1[ [Schwab et al.||2015 l. These formation mechanisms provide the simplest 

explanation for both the observed short spin periods P to as small as P = 1.39 ms, and 
superstrong magnetic fields with Bq ~ 10'^ G. 

Most radio pulsars are solitary. Of the more than 2400 pulsars known by the end of 2014, 
only about 230 were members of binary systems^Even in binary systems, mass transfer 
from the companion star to the neutron star is negligible. The radio luminosities of pulsars 
are low relative to the sensitivities of even the largest radiotelescopes, so that our catalog 
of pulsars is not complete even to a distance of a kpc. Because the Milky Way is an order 
of magnitude larger, we can observe only a small fraction of “active” pulsars. Because the 
duration of the active life of pulsars is small, the total number of extinguished pulsars in our 
Galaxy must be about 10*-10^ (Manchester et al. 2005 l 


2 Theoretical Overview 


2.1 Early Pulsar Paradigm — Vacuum Dipole 


The basic physical processes determining the observed activity of radio pulsars were un¬ 
derstood almost immediately after their discovery ( |Pacini|1967[|Gold|1968[ l. In particular, 
it quickly became clear that the highly-regular pulsed radio emission that gives rise to their 
name is related to the rotation of neutron stars. Furthermore, it was evident that radio pulsars 
are powered by the rotational energy of the neutron star, and the mechanism of energy re¬ 
lease is related to their superstrong magnetic fields, with Bq ~ 10*^ G. The Larmor formula 
for energy loss form a magnetic dipole provides an estimate of energy losses ( [Landau &| 
|Lifshitz|1989^ : 


Wtot = —ItQQ 

1 BlQ^R^ . 

~ 7-5— sm^^ 


(1) 

(2) 


where 4 ~ MR^ is the moment of inertia of the neutron star, x is the angle between the 
magnetic dipole axis and the spin axis, and Q = iTijP is the angular velocity of neutron star 
rotation. Finally, the strength of the magnetic field at the polar cap is Bq. 

For most pulsars, energy losses range from lO^’-lO^'* erg s“* and can reach 10^*^-10^® erg s 
for very young, fast pulsars, such as the Crab and Vela pulsars. These energy losses corre¬ 
spond to the observed spin-down rate dPjdt ~ 10“*^, or to the spin-down time tq = P/2P ~ 
1-10 Myr. 

' See ATNF catalog: http://www.atnf.csiro.au/people/pulsar/psrcat/ 
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After the measurement of the rotational slow-down P of the Crab pulsar ( |Richards~& 
|Comella|1969^ , it was quickly realized that: 

• the rate of the energy loss of the rotating neutron star Wtot ~ 5 x 10^* erg s coincides 

with the power required to illuminate the Crab Nebula ( |Gold|1969| l, and 

• the dynamical age of the radio pulsar td = P/2P ~ 1000 years coincides with the explo¬ 

sion of the historical supernova AD 1054 that brought the Crab Nebula into existence 
( [Comella et al.|1969^ . 

These associations cemented the identification of pulsars as rotating neutron stars. In con¬ 
trast to these phenomena, radio emission amounts to only 10 “^- 10 “^ of total energy losses. 
For most pulsars this corresponds to 10^®-10^® erg s“*, 5-7 orders of magnitude less than 
the luminosity of the Sun. 


2.2 Electron-positron generation 


IGoldreich & Julia^ ( |1969[ l showed shortly after the discovery of pulsars that a pulsar’s ro¬ 
tating magnetic field will acquire a corotating charge density that opposes induced electric 
fields and J x B forces. As |Sturrock| ( [1971| l quickly realized, individual photons can gen¬ 
erate electron-positron pairs when they cross lines of the magnetic field, by the process 
7 -H —> e'^ + e~ + B. The photon energy must exceed the threshold The proba¬ 

bility per-unit-length for conversion of a photon with energy By far above this threshold 
propagating at an angle of 6 to the magnetic field B is (|Berestetsky et al.|1982^ 


3 Vs / 8 Br, mgC^\ 

I 6 V 2 3Bsin0 By ) 


(3) 


Here, the characteristic value Bn = jeti ~ 4.4x 10*^ G is the magnetic field for which the 
energy gap between two Landau levels reaches the rest energy of an electron: Tiub = m^c^. 
As gamma-quanta are radiated by particles moving along the curved magnetic field lines, 
one can evaluate the photon free path as ([Ruderman & Sutherland] 1975|l 



Bn 

~B ~B^' 


(4) 


Here is the curvature radius and A s; 20 is the logarithmic factor. As, ly <K R for high 
enough photon energy, the vacuum magnetosphere of a neutron star with magnetic field 
Bq ~ 10*^ G is unstable to the generation of charged particles. 

In the very strong magnetic field of the neutron star, charged particles can move only 
along magnetic field lines. Therefore two substantially different regions must develop in 
the pulsar magnetosphere: regions of open and closed magnetic field lines (see Figure [T] 
1^. Closed magnetic field lines do not intersect the light cylinder, where co-rotation speed 
equals that of light, at radius Rl = cjQ (~ 10*° cm for ordinary pulsars). Particles on these 
field lines turn out to be captured. Open field lines intersect the light cylinder, and parti¬ 
cles on these field lines can travel to infinity. Consequently, plasma must be continuously 
regenerated near the magnetic poles of a neutron star (see Figure[^. 

In addition to the primary plasma generated by individual photons and the magnetic 
field, as discussed above, a secondary plasma forms from the longitudinal electric field 
(which accelerates particles up to energies high enough to radiate hard y-quanta), as first 
indicated by |Sturrock| (|1971|l and then studied in more detail by [Ruderman & Sutherland] 
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Fig. 1 Vacuum dipole model of a pulsar, 
showing the region of closed magnetic field 
lines and the light cylinder. The angle be¬ 
tween the rotation and magnetic poles is 
often called the “inclination angle”; and the 
minimum angle between the line of sight 
and the magnetic axis is sometimes called 
the ’’impact angle”. Sometimes x is denoted 
by a, and { by J3. See|Lyne & Manchester! 

(TW8) . 


Fig. 2 Intersection of the cone of the last open field 
lines with the surface of the neutron star, showing the 
polar cap. The vector m is the axis of the magnetic 
dipole field, and Ohs points toward the observer. 


( |1975| l, as well as by Eidman’s group ( [AEber et al.|1975y The continuous escape of parti¬ 
cles along the open field lines leads to formation of a strong electric field along the mag¬ 
netic field. This longitudinal electric field forms in the vicinity of the magnetic poles. The 
secondary plasma-generation condition determines its height. Another model, based on the 
assumption of free particle escape from the neutron star surface,was first studied by Arons’ 


group <|Fawley et aL|1977[ Scharlemann et aL| 

1978| Arons & Scharlemann|1979|), and re- 

cently in more detail by Istomin & Sobyanin 

(|2011a|bi; Medin & Lai (2010i; Timokhin| 

l|2010|l;|Chen & Beloborodov i 20131 and|Timokhin & Arons 1 2013b. 


2.3 Hollow-cone model 

The hollow-cone model ( [Radhakrishnan & Cooke|1969[|Dyks et al.|2004| l explains the basic 
observed properties of radio emission in the context of the above particle generation pro¬ 
cesses, without reference to a microphysical model for that emission. This model, already 
proposed at the end of the 1960s, perfectly accounts for the basic geometric properties of the 
radio emission. This model proposes that outflowing plasma launches radio emission tan¬ 
gent to open magnetic field lines at a particular altitude above the surface of the neutron star. 
The characteristic frequency of radiation may depend on altitude: the “radius-to-frequency 
mapping” l |Ruderman & Sutherland] 1975 [ I . Plasma density and geometry of open field lines 
define a “directivity pattern”. The observed average pulse is a cut across this directivity 
pattern. 
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Fig. 3 Structure of the particle generation 
region. The primary particles are acceler¬ 
ated along the curved magnetic field lines 
and begin to radiate hard y-rays. These cur¬ 
vature photons {dotted lines) propagating in 
the curved magnetic field reach the parti¬ 
cle generation threshold and create electron- 
positron pairs. Taken from|BeskinH199^. 


Fig. 4 The hollow cone model. If the intensity of the 
radio emission is directly connected with the outflow¬ 
ing plasma density, radio emission must decrease near 
the magnetic pole. Consequently, we expect a double 
profile when the line of sight passes near the magnetic 
pole, and a single profile when it passes further away. 
Taken from |Beslanl)1999) . 


Secondary particle generation is impossible in a nearly rectilinear magnetic field be¬ 
cause, first, little curvature radiation is emitted; and second, photons emitted by relativistic 
particles propagate at small angles to the magnetic field. Therefore, as shown in Figure]^ in 
the central region of the open magnetic field lines, a decrease in secondary plasma density 
is expected. If we make the rather reasonable assumption that radio emission is less when 
the outflowing plasma density is less, the intensity of radio emission must decrease in the 
center of the region of open field lines, corresponding to the center of the directivity pattern. 
Therefore, if without going into details]^ we should expect a single (one-hump) mean pro¬ 
file in pulsars in which the line of sight intersects the directivity pattern far from its center 
and the double (two-hump) profile for the central passage. This is precisely as observed in 
reality ( |Lyne & Graham-Smith|1998l l. 


3 Observational Overview 


Pulsars take their name from their remarkably stable periodic emission. The rotation fre¬ 
quency of the pulse train is the angular velocity of the neutron star. Folding the observed 
pulse train at this fundamental frequency yields an average pulse profile. In most cases this 
profile is extremely stable, both in form, and in arrival phase at the rotational frequency. 
This stability allows for precision timing of pulsars, with remarkable applications in struc¬ 
ture and evolution of stellar systems containing pulsars, and in tests of special and general 
relativity ( |Camenzind|2007[ l. The stability of the mean profile suggests that rotation carries 
the line of sight through a beam of emitted radiation locked to the surface of the neutron star. 


^ Actually, the mean profiles have a rather complex structure, see e.g. 
|Graham-Smith||l998^ 


Rankin 

1983 

1990 ; 

Lyne & 



























6 


V.S. Beskin et al. 


and that relatively permanent features of the neutron star and its co-rotating magnetosphere 
determine the shape of that beam. 

However, pulsar emission shows a remarkable degree of variability on all timescales, 
extending from nanoseconds to months or years. The stable pulse profiles that characterize 
that stability appear only after 100 or more pulses are added together, for pulsars strong 
enough to detect variability of emission. Indeed, [Popov et aL] ( |2006| l suggest that individual 
micropulses are the “atoms” of pulsar emission. For those who seek to understand emission 
processes of pulsars, as well as those who merely wish to exploit pulse stability for other 
scientific goals, pulse variability can provide crucial insights. 

Observations of pulsars occupy a multi-dimensional space. The fundamental observ¬ 
ables include intensity and polarization of electromagnetic radiation, as function of time in 
pulse phase and over many pulses. These fundamental observables show both determinis¬ 
tic and random properties, with random properties in particular showing variations over all 
timescales. At radio wavelengths, pulsar spectra are nearly power-law, but comparisons of 
pulse shape and structure among wavelength ranges yields great insight into emission ge¬ 
ometry and processes (see, for example, [Shearer et al.|2003) [Lommen et al.|2007{ [ Harding[ 
[et al.|2008l|Strader et al.|2013^ . 

Among the important quantities derived from observations of radio emission from pul¬ 
sars are the spindown rate, polarization as a function of pulse phase, size of emission region, 
and evolution of angle between the spin and magnetic dipole. 

4 Magnetosphere of an Oblique Magnetic Rotator 

A pulsar represents an elegant problem in electrodynamics: a rotating, conducting sphere 
with a dipole magnetic field ( [Beskin & Zheltoukhov|20l4[ l. This simple picture is compli¬ 
cated by the necessity of a corotating charge distribution, the roles of open and closed field 
lines, and energy transport by an outflowing wind ( [Goldreich & Julian|1969^ . As we sum¬ 
marize in this section, models have progressed from analytic studies of aligned rotators with 
simple magnetic field configuration and massless charges to self-consistent models includ¬ 
ing oblique magnetic fields, realistic particle masses, and a range of length scales. 

4.1 Current Losses 

If the pair creation process is sufficiently effective, magnetic dipole radiation will not carry 
energy away from the rotating neutron star, because the plasma that fills the magnetosphere 
fully screens any low-frequency radiation from the neutron star ( [Beskin et al.|198^|1993 [ 
[Mestel et al.[1999[ l. However, in this case, electric currents extract rotational energy from the 
neutron star, through the Ampere force of one current on another. The currents in question 
are those along magnetic field lines in the magnetosphere and across the pulsar’s polar cap, 
acting together with those responsible for the magnetic field of the neutron star. Just as in 
the case of magnetic dipole radiation, energy release from the rotating neutron star is related 
to the electromagnetic energy flux given by the Poynting vector, and the total energy losses 
can be again estimated using the Larmor formula, ([^. 

The braking torque K of the Ampere force results in the following time evolution of the 
angular velocity Q and the inclination angle x- 


IjQ = cosx + sin;t^, 
IrOx = Ki. cosx - K\\ sinx. 


(5) 

( 6 ) 
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Fig. 5 Schematic view of the axisymmetric polar cap 
showing magnetospheric current density (open ar¬ 
rows), surface currents. Ampere force on surface cur¬ 
rents, and braking torque. Here only the symmetric 
current is present. Taken from|BeskinH199^. 



Fig. 6 For the orthogonal rotator only antisymmetric 
current 4 (i-e-. the current having different direction 
in the north and south parts of the polar cap) takes 
place. The structure of the surface currents within 
the polar cap and along the separatrix is also shown. 
Taken from |Beskin et al.|^2013) . 


where two components of the torque K parallel and perpendicular to the magnetic dipole m 
can be written in the form l |Beskin et al.|1993l t 


^11 = 

K, = 


-C||- 




-Cx- 


BlQ^R^ 


QR 

c 


(7) 

( 8 ) 


Here the coefficients cy and Cj_ are factors of the order of unity dependent on the profile of 
the longitudinal current and the form of the polar cap. 

The scalar current from the polar cap i has been divided into symmetric and antisym¬ 
metric contributions, 4 and 4, depending upon whether the direction of the current is the 
same in the north and south parts of the polar cap, or opposite. For an axisymmetric rotat¬ 
ing neutron star (x = Q, Figure [^, we have 4 = 0 and 4 = 1 ( [Goldreich & Julian|1969^ . 
Conversely, for the orthogonal rotator (Figure]^ we have 4 = 0 and 4 = 1- Here we apply 
normalization to the Goldreich-Julian current, /qj = nR^ja, where Rq ~ RiQRjcy^^ is the 
polar cap radius, and ^'gj =< |n • B| > jin (with scalar product) is the mean current density 
within the polar cap. Note that for 4 ~ 4 ~ 0 0 imply that: 


^x 



(9) 


Therefore, A'x K\\. We will use these expressions in the following sections. 

If we suppose that, in reality, the longitudinal current j is determined by the local charge 
density pcj = -£l ■ B/2;rc, and note that pcj is proportional to cos;t^ in the vicinity of the 
polar cap, one can write down 


4 = cos^, 

4 = it sin^. 


( 10 ) 

( 11 ) 
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Fig. 7 The Michel split-monopole solution, 
in which electric field Eg has magnitude 
equal to the toroidal magnetic field This 
solution requires a conducting current sheet, 
to close electric currents (outlined arrows). 

Taken from |Beskin| l | 199^ . 

Consequently, the relations can be rewritten in the form 

IrQ = K^ + lK^-K^]sm^X, (12) 

hQx = - K^] sin;t-cos^. (13) 

As we see, both expressions contain the factor [K^ - This implies that the sign of^ 
is given by the ;^^-dependence of the energy losses ( [Philippov et al.|2014[ l. In other words, 
the inclination angle x will evolve to 90° (to counter-alignment) if the total energy losses 
decrease for larger inclinational angles, and to co-alignment if they increase with inclination 
angle. 

Because the plasma filling the pulsar magnetosphere is secondary (in other words, it is 
produced by the primary particles accelerated by the longitudinal electric field), at any point 
out to the light cylinder, the energy density of the electromagnetic field must be much larger 
than the energy density of the magnetospheric plasma. For the same reason, energy transport 
is given by the Poynting vector (see Figure [^. 


Con- 


Fig. 8 Magnetospheric structure in the model of 
|topouIos et al.|{T9991 . The last open field line is as¬ 
sumed to coincide with the equator. 




4.2 Split-Monopole Model 

The remarkable analytical solution found by |MicheT] ( |1973[ l serves to illustrate the transport 
of energy by Poynting flux. In the force-free approximation, when massless charged particles 
move radially with the velocity of light, and with the Goldreich-Julian current density jci = 
Pg]C, a split-monopole magnetic field is the exact solution to the Maxwell’s equations, both 
inside and beyond the light cylinder (see Figure [^. The monopolar magnetic field is split 
so that the magnetic flux converges in the southern hemisphere and diverges in the northern 
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one. In this solution. Ampere forces from longitudinal currents along magnetic field lines, 
and from corotation currents from rotating charge density, are fully compensated. 

In the Michel split-monopole solution, the electric field E has only a 0-component, and 
is equal in magnitude to the toroidal component of the magnetic field: 

(QR\ R 

5^ = Efl =-Sol —1 - sine. (14) 


At distances larger than the light cylinder radius, this magnetic field becomes larger than the 
poloidal magnetic field Sp = So(S/r)^. On the other hand, in this solution the total magnetic 
field remains larger than the electric field everywhere, so that the fields form electromagnetic 
waves only at infinity. 

Because magnetic flux converges in the lower hemisphere and diverges in the upper one 
in the split-monopole solution, a current sheet must lie in the equatorial plane (see Figures]^ 
[^. This sheet closes the longitudinal electric currents elsewhere in the magnetosphere. This 
structure of the magnetic field and current sheet has been confirmed numerically (|Contopou 


|2006[ |McKinney [2006 


los et: 


alJ1999 Ogura & Kojima|2003t|Gruzinov|2005][Ki 


omissarov 


Timokhin|20061. 


Using Eqn. ( |14[ l, one easily finds that the Poynting vector S = (c/47r)E x B is: 


(15) 


This implies that the energy flux is concentrated near the equatorial plane. This 0-dependence 
of the energy flux is used by many authors (|Bogovaiov & Khangoulyan |2002[ Komissar^ 
& Lyubarsky 2003[ l. On the other hand, at large distances r » R]^, Ingraham 11973 I and 


Michel 


11974 I found another asymptotically radial solution, with Eg{6) = B^p{6), resulting in 


a radial Poynting vector with arbitrary 0-dependence. In all of these solutions, the relation 

(16) 


s{e)ozBl{e)s.m^e 


is valid. 

|Bogovalo^ ( | 1 999[ l generalized the split-monopole model, showing that in the force-free 
approximation the “inclined split monopole field” is a solution of the problem as well. In 
this solution, 

IQR\ R 

= Efl =-So II - sm 0 sign© (17) 

and Sp = Bo(Rlr)^ sign©, where 

© = sin^ sin 0 sin(^ — Qt + Qrjc) + co&x cos 0. (18) 


In this case, within the cones 0 < njl - x, n - 6 < njl — x around the rotation axes, the 
electromagnetic field is not time dependent; whereas in the equatorial region, the electro¬ 
magnetic fields change the sign at the instant © = 0. In other words, the condition © = 0 
defines the location of the current sheet. We stress that the expression for the shape of 
the current sheet remains true for the other radial asymptotic solutions, with EgiG) = S^(0) 
but arbitrary 0-dependence, as well l |Arzamasskiy et al.|2015a| l. Numerical simulations ob¬ 
tained recently for the oblique force-free rotator confirm this conclusion as well (Spitkovsky 


2006 Kalapotharakos & Contopoulos(200^ [Kalapotharakos et al.|2012[|Tchekhovskoy 

al.|2013 Philippo'v et al.|20l4 1. 
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4.3 Magnetohydrodynamic Models 

As was already stressed, recently numerical simulations have become possible that can sim¬ 
ulate the structure of plasma-filled magnetospheres from first principles. |Contopoulos et 31] 
( |1999[ l found an iterative way to do this and obtained the first solution for an aligned force- 
free pulsar magnetosphere that extended out to infinity (see Figure [^. Their results were 
subsequently verified by other groups within force-free and magnetohydrodynamic (MHD) 


vomrssarovl. 


Parfrey et al.|2012 

[Ruiz et al. 2014 1 as well as using particle-in-cell (PIC) approach ( Philip-j 

pov & Spitkovsky 

2014 Chen & Beloborodov|2014[|Cerutti et al.|2014[|Belyaev|20I4 1 . 


Using the force-free approximation, he found that pulsar spindown luminosity increases 
with increasing obliquity angle, which is the angle between the rotational and magnetic 
axes. The spindown obtained in such force-free and MHD models is well-described by 


~ fl^aligned(l "k siu ^), 


(19) 


where VTaiigned = nP'Q* I is the spindown luminosity of an aligned plasma-filled pulsar 
magnetosphere, and m = BoF? j2 is the magnetic dipole moment of the pulsar. More re¬ 


cently, these results were confirmed using time-dependent 3D force-free i Kalapotharakos 
& Contopoulos|2009 Petri|201^ [Kalapotharakos et alj20T^, MHD I jTchekhovskoy et^ 


2013|l, and PIC (Philippov et al.|2014i studies. Figureshows a vertical slice through the 


results of a 3D MHD simulation of an oblique pulsar magnetosphere with obliquity angle 
X = 60°. One can clearly see the closed zone that extends out to the light cylinder located 
at |x| = beyond which starts a warped magnetospheric current sheet, across which all 
field components undergo a jump. The structure of this current sheet is presently poorly un¬ 
derstood, in particular it is not known if in the perfect conductivity limit the magnitude of 
the magnetic field in the sheet vanishes. 

What causes this increase of spindown luminosity with the increase of pulsar obliquity? 
It turns out that there are two factors: (i) an increase in the amount of open magnetic flux, 
which accounts for about 40% of the increase, and (ii) redistribution of open magnetic flux 
toward the equatorial plane of the pulsar magnetosphere, which accounts for the remaining 
60%. 

In fact, the spindown trend ( |19| l can be reproduced via a simple toy model. Suppose that 
the magnetic field that reaches the light cylinder in an oblique rotator, with inclination angle 
X, retains the dipolar structure at r = ro » i?LC> 


Br = —^singm, (20) 

yL 

where 6^ is the magnetic colatitude, or the angle away from the magnetic axis, 

0 ni = arccos(sin^sin0cos^-H cos6cos;^f). (21) 

How would the pulsar spindown change if we kept the total open magnetic flux, (Popen = 
nr^Bo, fixed, and inclined the pulsar, i.e., increased^? To find this out, let us first compute 
the angular distribution of ^-averaged B^: 


(Bj)y; = bI{0.5 sin^;(f sin^ 6 + cos^ 6cos^;y). 


(22) 
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Fig. 9 Slice through the m-H plane of a relativistic MHD simulation of an oblique pulsar magnetosphere 
(obliquity angles = 60°) taken after 3 rotations. Solid lines show field lines as traced in the image plane. 
Into-the-plane magnetic field component, is shown with color (red — into plane, blue - out of plane). 
Taken from |Tchekhovskoy et al.| l [2013^ . 


Now, making use of the fact that ^ Eq = -BfOrsmOjc and the radial Poynting flux is 
Sr = cEgB^jAn ^ {BrQr^in&f-jAnc, we obtain ( Tchekhovskoy et al.|2015 i: 


rr 0 A 

'^loyix) = jj S,.doj= + sin^/T) 


(23) 


where the integral is over, e.g., a sphere of radius tq. Clearly, if the total magnetic flux tPopen 
is held constant, the nonuniformity in the surface distribution of magnetic flux causes an 
enhancement in spindown losses at higher inclination angles, consistent with the numerical 
simulations (see Eqn. [T^ . In the simulations, we find that the magnetic flux itself is an 
increasing function ofx, t^open oc (1 + 0.2 sin^;^')|^so the two effects - of the non-uniformity 
of the open magnetic flux and the change in the amount of open magnetic flux - have very 
similar inclination-dependences. 

In reality, the structure of the plasma-filled magnetosphere is of course more complex 
than given by Eqn. ( |20[ l, but the qualitative effect is the same: the inclination of the magnetic 
axis relative to the rotational axis leads to the shift of the peak of |Srl away from the axis 
and toward the equatorial plane and an increase in the spindown luminosity (|Tchekhovskoy| 
let al.|2015| l. 


More recently, PIC models have been developed and, in those cases when the magne- 
tospheric polar cascade is efficiently operating and is able to fill the magnetosphere with 
abundant plasma, are in agreement in the amount of spindown and large-scale dissipation 


as in MHD simulations (see e.g. 

Philippov & Spitkovsky 

2014HChen & Beloborodov|2014| 

ICerutti et al.|2014||Belyaev|2014| Philippov et al.|2014 

1 . Interestingly, if a mechanism of 


3 


This relation coincides exactly with one obtained by 


Beskin et al. 


1993 


analytically. 
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pair formation operates only near the surface of the star, aligned pulsar magnetospheres 
in PIC simulations do not reach a force-free state ( |Chen & Beloborodov||201^ . In fact, 
PIC simulations, into which simplified physics of the polar cascade was included, show 
the development of the polar cascade and of a force-free-like magnetosphere only for high 
inclinations, a > 40° ([Philippov et al.|2014^. 


5 Observations: Energy Loss from Pulsars 

5.1 Spindown 

In principle, the time rate of change of pulse period P is easy to measure. Because indi¬ 
vidual pulses can be numbered, period and period derivative are among the fundamental 
parameters of a timing model. Period derivative is easily associated with the loss of rota¬ 
tional kinetic energy via electromagnetic radiation and particle wind. The Larmor formula 
for magnetic dipole radiation then directly associates energy loss with the magnetic moment 
of the neutron star. This provides a characteristic scale. 

Pulsars with periods longer than a fraction of a second show timing noise: random vari¬ 
ations of pulse arrival time that change slowly with time ( jHelfand et al.|19 ^. These vari¬ 
ations are most extreme for the young Crab and Vela pulsars ( Boynton et al.|[T972[ |Lyne| 
|& Graham-Smith| 1998[ [Scott et al. [2003 [[Dodson et al.|2007| l. Among millisecond pulsars, 
B1937-1-214 shows timing noise, but other millisecond pulsars may not ([Kaspi et al.|1994| 
[Cognard et al.|19^ . 

Several pulsars show clear variations in spindown rate associated with changes in pulse 
properties. The radio pulsars B1931 h-24, I1832-h 0029, and J1841-0500 intermittently switch 
between an “on” radio-loud state in which they appear as ordinary radio pulsars, and an “off” 
state in which no radio emission is detected. The spin-down rate is higher in the “on” state 
than the “off” state, by a factor of /on^oft = ^on /Pofi ~ 1.5 for B1931-H24 (|Kramer et al.| 
[2006[ l and J1832-H0029 ( [Lyne[2009] l, and/on->off ~ 2.5 for J1841-0500 ( [Camilo et al.[2012] ir 
The gamma-ray pulsar J2021-t4026 [Allafort et al.[ ( [2013[ l displays two states with intensities 
different by 20% and with distinct pulse profiles, each associated with a different spindown 
rate: fon->ofi = 1.04. Pulsar B0919-H06 shows quasiperiodic variations between two states 
with different spindown rates and different pulse profiles ( [Perera et al.[[2015| l. [Lyne et aT] 
l [2010[ l propose that the phenomenon of intermittency is quite general: they find that timing 
noise for six pulsars can be expressed as the superposition of two states, characterized by 
distinct pulse profiles and spindown rates, with rather rapid changes between states. From 
these discussions it is clear that magnetospheric structure affects spindown, as one would 
suspect from theoretical considerations discussed in Section 2.3 above. 

[Kramer et al.[ ( [2006[ ( (see also [Beskin & Nokhrin^[2007^ proposed that two distinct 
magnetospheric states lead to the observed difference in spindown rates. They associated the 
“off” state with a magnetosphere depleted of charge, and the “on” state with magnetospheric 
currents sufficient to produce the observed change in spindown. [Li et al.] (2012 1 observe that 
the simplest model for the “on” state is the force-free magnetosphere i Spitkovsl^[2006[ l, 
which exhibits spindown rates at least three times that of a vacuum dipole. They suggest a 
modified picture where the “on” state is the force-free magnetosphere, and the “off” state 
has no charge on open field lines, but carries the Goldreich-Julian charge on closed field 
lines. This leads to ratios /on^otf = 1.2 to 2.9 for inclination angles of;^^ > 30°. Smaller 
inclinations lead to larger fon->oS- 
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5.2 cr-Problem 


Thus, we see that all analytical and numerical force-free models of the pulsar magnetosphere 
demonstrate the existence of an almost-radial highly-magnetized wind, flowing outward 
from the pulsar magnetosphere. On the other hand, observations show that most energy far 
from the neutron star must be carried by relativistic particles ( [Kennel & Coroniti|1984a|b) l. 
For example, the analysis of the emission from the Crab Nebula in the shock region located 
at a distance of ~ 10*^ cm from the pulsar in the region of interaction of the pulsar wind 
with the supernova remnant definitely shows that the total flux V 14 m of the electromagnetic 
energy in this region is no more than ~ 10“^ of the particle energy flux ITpan. Thus, in the 
asymptotically far region of pulsar models, the Poynting flux must be completely converted 
into an outgoing particle flux before reaching the reverse shock at distances of ~ 0.1 pc. 
Axisymmetric numerical models of jets from radio pulsars are constructed exactly under 
this assumption ( [Kirk et al.|2009[ and references therein). 

The transformation from Poynting flux to particles apparently occurs much closer to the 
neutron star, at distances comparable to the size of the light cylinder. This is evidenced by 
the detection of variable optical emission from companions in some close binary systems 
involving radio pulsars (Fruchter et al.||1988[ Kiflkarnirt^alJ|]^88|^ Fruchter et al. | 


Ryba & Taylor|199T| Stappers et al.j 19961 Roberts 2011 Pallarica et al.|2012y Romani et al. 


2012[ ^plan et al.|2013HBreton et al.|2013 I. This variable optical emission with a period 

equal exactly to the orbital period of the binary can be naturally related to the heating of the 
companion’s surface facing the radio pulsar. It was found that the energy reradiated by the 
companion star almost matches the total energy emitted by the radio pulsar into the corre¬ 
sponding solid angle. Clearly, this fact cannot be understood either in the magnetic-dipole 
radiation model or by assuming a Poynting-dominated strongly-magnetized outflow, since 
the transformation coefficient of a low-frequency electromagnetic wave cannot be close to 
unity. Only if a significant fraction of the energy is related to the relativistic particle flux 
can the heating of the star’s surface be effective enough. Moreover, eclipses of the double¬ 


pulsar system show effects of the particle wind from one object impinging upon (Lyne et al. 


|2004[penet & Ransom|20(34)|McLaughlin et al.|2004[|Lyutikov|2004[[Demorest et al.|2004^ 

Therefore, the so-called cr-problem — the question as to how the energy can be converted 
from electromagnetic fields to particles in the pulsar wind — remains one of great unsolved 
problems of modern astrophysics. We note that the cr-problem appears to be rather general 
and in addition to neutron-star powered outflows it applies to black-hole powered, collimated 
outflows known as astrophysical jets, such as in the magnetically-arrested disk (MAD) sce¬ 
nario ( Tchekhovskoy et ak] 20]TJ^ Tchekhovskoy^^cKmneyjj20]^ jZamaninasab et aL] 
|2014[[Mziarski et al.|2014[|(jhisellini et al.|2014[|Tchekhovskoy|2015| . Theoretical models 
suggest that the jets accelerate roughly up to the equipartition between the magnetic and 
kinetic energies, beyond which the acceleration slows down dramatically, locking in a sub¬ 
stantial fraction of energy in the magnetic form (Tchekhovskoy et al.|2009[ [Komissarov etj 
[al.|2009t [Lyubarsky|20 1 0[ however, see [Tchekhovskoy et al.|20i0 1 . 


6 Theory: Polarization and Refraction of Radio Emission 

6.1 Polarization 

Pulsar emission is usually highly linearly polarized, with a small fraction of circular polar¬ 
ization. Like the mean profile, the profiles in polarization states are stable and are charac- 
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teristic of the pulsar. This long-term stability of the mean properties indicates that the pulse 
arises as a cut through a radiation cone, with properties that are set by stable properties of 
the underlying neutron star. It is widely assumed that the polarization is determined by mag¬ 
netic fields in or above the emission region. Those magnetic fields, in turn, are anchored in 
the solid crust of the neutron-star ( |Manchester| 1995 | l . 

Rapid swings of the position angle of linear polarization through the pulse, first observed 
in the Vela pulsar by |Radhakrishnan et al.| ( [T96^ , suggest that a vector fixed in the frame of 
the rotating star influences the direction of linear polarization, a geometric inference known 
as the rotating-vector model. [Radhakrishnan & Cooke| ( |1969| l proposed that this vector is 
the magnetic pole of the pulsar’s nearly-dipolar magnetic field; this physical interpretation 
is known as the magnetic-pole model. 


Radiotelescopes can measure polarization properties of individual pulses for a number 
of strong pulsars. Such studies indicate the presence of orthogonal modes, with polarization 


differing by 90°, and intensities varying from pulse to pulse (Manchester et aL|1975 

|B acker 

et al.|1976HCordes et al.|1978||Backer & Rankin] 1980||Stinebring et al.|1984a|b[ 

McKin- 

non|2003|l. For most pulsars, present radiotelescopes can determine only average polariza- 


tion properties; nevertheless the presence of two competing orthogonal modes can explain 
the observed departures from the characteristic pattern, for most of these weaker pulsars. 
Thus, the rotating-vector model, with two orthogonal linearly-polarized modes, success¬ 
fully describes the characteristic swing of the angle of polarization with pulse phase for 
most pulsars, across a wide range of pulsar parameters and despite observational selec¬ 
tion effects l |Rankin|ri983| |1986[ |Lyne & Manchest^|1988| l. The two modes are usually 
interpreted as the X-mode, with wave electric field perpendicular to stellar magnetic field 
{Ew ± 5o); and the 0-mode, with a component of electric field parallel to stellar field 
{Ev/\\Bq). Both modes appear to be present, at some level, for all radio pulsars. 

Work to model the polarization properties of pulsars in more detail, including the circular- 
polarized profile, have led to mapping of the polarization properties on the Poincare sphere 
describing the Stokes parameters ( |McKinnon|2009[[Chung & Melatos|201 la|b| l. These show 
a rich variety of patterns, with greater modulation of polarization being indicative of more 
complex patterns. Analysis of these patterns suggest emission, or refractive scattering within 
the pulsar’s light cylinder. For some pulsars, the emission, or reprocessing region is inferred 
to lie at altitudes of 10 to 40% of the light-cylinder radius. 


6.2 Rotating Vector Model 


The standard relation for the rotating vector model describes variation of the position angle 
of polarization tp in the mean profile, under the assumption that the the hollow-cone model 
is valid. In other words, it assumes that all absorption is absent, and that the magnetic field is 
dipolar in the emission region, precisely where the polarization is determined. This relation 
takes the form: 


\p = arctan 


sin;(f sin<^ 

sin;t' cos f cos - sin ^ cos;^^ 


(24) 


Here, once again, is the inclination angle of the magnetic dipole to the rotation axis, f is the 
angle between the rotation axis and the direction toward the observer, and ^ is the phase of 
the pulse. Figure[T^illustrates the geometry. Equation \2A) has been used for many years in 
estimating the pulsar inclination angle, which is a very important parameter for determining 
the structure of the magnetosphere. Aberration and retardation effects ( Blaskiewicz et al.| 
|1991|l have been included in only some studies (|Mitra & Li|2004[|Krzeszowski et al.|2009^. 
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Fig. 10 Geometry for the rotating vector model. The emission surface is shown as a golden band. As the 
pulsar rotates, the line of sight to the observer traverses a circle of constant latitude, producing an observed 
pulse each time the emission surface crosses the line of sight. The dipolar magnetic field, projected onto 
the radius of the emission surface, radiates from the magnetic pole m. The instantaneous polarization of the 
observed radiation is parallel to the field lines, at position angle if/ at the observer. 


The rotating vector model, extended further with the “hollow cone” model, is based 
on the following three assumptions (see, e.g., [Manchester & Taylor| ( |1977[ l): the formation 
of polarization occurs at the point of emission; radio waves propagate along straight lines; 
and cyclotron absorption can be neglected. But all these assumptions turn out to be incor¬ 
rect. [BamMj^^o^ (|T98^ showed that in the innermost regions of the magnetosphere, 
the refraction of one of the normal modes is significant. After publication of the work of 
Mikhailovskii’s group ( [Mikhailovskii et al.|1982| l, it became clear that cyclotron absorption 
can significantly affect the radio emission intensity. The influence of the magnetosphere 
plasma on variation of the polarization of radio emission propagating in the internal regions 
of the magnetosphere also must not be neglected ( [Petrova & Lyubarskii|2000l l. 

The “limiting polarization” is the most important effect of magnetospheric propaga¬ 
tion. Radio emission in the region of dense plasma consists of a superposition of normal 
modes: in particular, the principal axes of the polarization ellipse must remain aligned with 
the magnetic-field direction in the picture plane. Polarization in the vacuum region is in¬ 
dependent of magnetic field. Hence, between the two lies a transition layer, past which the 
polarization is no longer affected by the magnetospheric plasma. For typical parameters of 
the pulsar magnetosphere, the formation of polarization occurs not at the emission point but 
at a distance of about O.IRl from it ( [Cheng & Ruderman|1979[[Barnard|1986[ l. Taking this 
effect into account should also explain the observed fraction of circular polarization of the 
order of (5-10)%. Therefore, a consistent theory of radio wave propagation in the magneto¬ 
sphere is required for a quantitative comparison of theoretical results on radio emission with 
observational data. 
















16 


V.S. Beskin et al. 


6.3 Propagation effects 

At present, the theory of radio wave propagation in the magnetosphere of a pulsar can be 


considered to provide the necessary precision (|Petrova|2006| |Wang et al.|2010t [Beskin &| 


|Philippov|201 1[ 2012 

Kravtsov & Orlov|1980[. Four normal modes exist in the magneto- 

sphere 1 

Beskin et al.| 

993|[Lyne & Graham-Smithl 1 998|l. Two of them are plasma modes 


and two are electromagnetic, which are capable of departing from the magnetosphere. An 
extraordinary wave (the X-mode) with the polarization perpendicular to the magnetic field 
in the picture plane propagates along a straight line, while an ordinary wave (the 0 -mode) 
undergoes refraction and deviates from the magnetic axis. An important point here is that for 
typical magnetosphere parameters, refraction occurs at distances up to 0 . 1 f?L, i-e., it can be 
considered separately from the cyclotron absorption and the limiting polarization. As shown 
in Figure [TT] the pulsar B0329-H54 shows both X- and 0-modes, with the 0-mode display¬ 
ing deviations from the rotating-vector model because of refraction ( [Edwards & Stappersj 
[2004l l. 

Based on the [Kravtsov & Orlov| ( [1980| l method, [Beskin & Philippov| ( [2012[ l have used 
such a theory of wave propagation in a realistic pulsar magnetosphere, taking corrections to 
the dipole magnetosphere into account (based on the results obtained by numerical simula¬ 
tion in |Spitkovsky[ ( [20()6) l), together with the drift of plasma particles in crossed electric and 
magnetic fields, and a realistic particle distribution function. The theory developed allows 
dealing with an arbitrary profile of the spatial plasma distribution, which may differ from 
the one in the hollow-cone model, because precisely the inhomogeneous plasma distribution 
leads to the characteristic ‘patchy’ directivity pattern ( [Rankin|1990[ l. 

The main result consists in the prediction of a correlation between the sign of the circular 
polarization (the Stokes parameter V) and the sign of the derivative of the change in the 
polarization of the position angle, if/, along the profile, where cj) is the phase of the 

radio pulse. For the ordinary mode, these signs must be opposite to each other, while for the 
extraordinary mode, they must coincide. FigurefTTjshows this pattern as well. As was noted, 
refraction of the ordinary wave leads to a deviation of beams from the rotation axis, and 
therefore the ordinary wave pattern should be broader than for the extraordinary wave. In the 
case of the ordinary mode, double radio emission profiles should mainly be observed, while 
single profiles should be observed in the case of the narrower extraordinary mode (jBeskinj 
|etal.[1993[ l. 

As was shown by [Andrianov & Beskin[ ( (20l0| , observations fully confirm the prediction 
of correlation between signs of V and d^/d<^. The analysis used over 70 pulsars with well- 
traced variation of the position angle and the sign of the circular polarization V, chosen from 
reviews of pulse profiles [Weltevrede & Johnston[ ( [2008| ; [Hankins & Rankin| ( [2010| l. Table[T] 
presents the results of the analysis. Pulsars with opposite signs of the derivative dtf/jdcp and 
the Stokes parameter V were placed in class O, while those with identical signs were placed 
in class X. As can be seen from the Table, most of the pulsars exhibiting a double-peaked 
(index D) profile indeed correspond to the ordinary wave, while most of the pulsars with 
single-peaked profiles (index S) correspond to the extraordinary wave. Moreover, the aver¬ 
age width of the radiation pattern for Od pulsars is indeed about two times larger than the 
average width of the radiation pattern for Xs pulsars. For the pulse width, the analysis used 
the width at the 50% intensity level VT 50 , normalized to the pulsar period P. The existence of 
a certain number of pulsars of classes Od and Xs should not give rise to surprise, because for 
central passage through the directivity pattern, independently of whether it corresponds to 
the 0-mode or to the X-mode, a double-peaked profile should be observed, while for lateral 
passage, a single-peaked profile should be observed. 
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Fig. 11 Polarization properties of the strong pulsar B0329+54. Upper panel: Average pulse profile, showing 
both total (black) and polarized (red) intensity plotted with pulse phase. Middle panel: Histogram of angle 
of linear polarization i}/ plotted with pulse phase. Dotted curves show polarization for the rotating-vector 
model. Red shows most common values, ranging through less-common green, blue, and violet, to white for 
not observed. One linear polarization (presumably X) tracks the rotating-vector model well, while the other 
(presumably O) shows large variations in polarization and timing. Lower panels: Polarization parameters on 
the Poincare sphere, for the 3 intervals of pulse phase indicated by red boxes at in the upper left panel. Each 
pair of disks shows 2 hemispheres. Linear polarization lies on the equator; pure right circular polarization is 
at the upper pole, pure left circular at the lower. This image from |Edwards & StappersH2004) is reproduced 
by courtesy of B. W. Stappers. 


Accurately taking propagation effects into account, [Andrianov & Beskin| t|2010[t; |Beskin| 
|& Philippov| l [20ir] t showed that such a variation of the position angle can be realized only 
under conditions of low plasma density or high mean particle energy. They found significant 
deviations from the standard relation of the rotating vector model (Eq. [24) 1 were obtained 
in the case of quite reasonable parameters that satisfy models for particle production: for 
example, a multiplicity ng/ncj ~ lO'^ and an average Lorentz factor y ~ 50. 
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Table 1 Statistics of pulsars with known circular polarization V and variation of position angle ijj. 


Polarization Mode O X 


Profile Type 

Single 

Double 

Single 

Double 

Class 

Os 

Od 

Xs 

Xd 

Number of Pulsars 

6 

23 

45 

6 

Normalized Pulse Width" 

6.8± 3.1 

10.7± 4.5 

6.5± 2.9 

5.3± 3.0 


a Normalized pulse width given as: VpVV'so deg) 


7 Observations: Polarization and Pulsar Size 


7.1 Pulsar Emission Region Size and Shift 


The radio emission regions of pulsars lie within the light cylinder, and so have angular 
sizes from Earth of nanoarcseconds or less. Resolving such an angle at radio wavelengths 
requires an instrument with an aperture approaching an AU, beyond the capabilities of even 
the longest VLBI baselines ( |Kardashev et al.|2013| l. However, radio-wave scattering by the 
dilute, turbulent interstellar plasma yields the required angular resolution and offer some of 
the information provided by a lens of that aperture. 

Interstellar scattering affects almost all astrophysical sources at decimeter wavelengths, 
and for many at shorter wavelengths. For most radio pulsar observations, scattering is “strong” 
in the sense that paths contributing to the electric field measured at the observer differ in 
length by many wavelengths. Hence, these paths behave like a corrupt lens ( |Gwinn et al.| 
|1998| l. The angular extent on the sky of these paths, 6 , delineates the “scattering disk”. The 
scale of variation of the impulse-response function at the observer, Siss, is the diffractive 
spot size of that aperture, ~ AjO. (Here, the subscript “ISS” indicates “interstellar scatter¬ 
ing”.) The scattering produces a random diffraction pattern in the observer plane with lateral 
scale 5 ISS - The effective resolution limit of the corrupt lens, at the source, is MS'iss, where 
the magnification factor M = DIR is equal to the distance D of the scattering material from 
the observer, divided by its distance R from the source. Interstellar scattering does not re¬ 
move information from the pulsar signal; rather, it adds a great deal of information about 
the paths taken. The challenge facing the observer is to extract the spatial information about 
the source from the scattered pulsar signal. 

Studies of the sizes of pulsar emission regions using interstellar scattering fall into two 
categories. One category relies upon the fact that if the emission point of the pulsar shifts 
across the pulse, the random image in the plane of the observer will undergo a reflex shift, 
as illustrated in Figure[T^ Proper motion causes a similar shift, but over time spans of many 
pulses. Correlation of the scintillation spectrum across pulse phases with later or earlier 
times yields the shift of the emission point ( Backer|1975[ Cordes et al.|1983t[Wblszczan &| 


|Cordes|1987[[Smirnova et al.|1996[ [Gupta et al.|1999[|Pen et al.|2014 1 . 

A second category invokes the decreased modulation for scintillation of an extended 
source. (“Stars twinkle, planets do not.”) The depth of modulation reveals the size of the 
I region (|Cohen et al.|1966t [Readhead & Hewish|1972[ [Hewish et al.|1974^ [G 


et al.||2012b[ [Johnson et al.|2012b l. More precisely, source size affects the distribution of 

flux density for a scintillating source. In strong scattering many different paths, with lengths 
differing by many radians of phase, contribute to the electric field measured at a point in the 
observer plane. The observer implicitly sums over these paths, so that the observed phase 
and amplitude have the character of a random walk. The optics of this effect are similar to 
those of the reflex shiff: different parts of the source produce shifted, incoherent diffraction 
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Fig. 12 Geometry for studies of pulsar structure by scattering. Radiation emitted by the pulsar at left travels 
to the scattering screen, where fluctuations in plasma density change its phase. The radiation arrives at the 
plane of the observer at right from along many paths, with different phases and amplitudes. Interference 
among paths produces a random diffraction pattern in the plane of the observer. When the source shifts to a 
different position, because of either proper motion or a rotational shift in the location of the pulsar emission 
region, the phases of the paths are modified, to produce a reflex motion of the diffraction pattern in the plane 
of the observer. 


patterns at the observer, who sums over them. Thus, finite source size affects the distribution 
of intensity at one antenna, or that of correlated flux density between the two antennas of an 
interferometer, principally by shifting the lowest and highest intensities toward the central 
part of the distribution ( |Scheuer|l968} [Gwinn|200T] [Johnson & Gwinn|2012a[|2013[ l. For 
realistic observations, the contributions of background noise, and of the noiselike statistics 
of the source itself, must be taken into account i jGwinn et al.|201 l||2012aHJohnson & Gwinn| 
|2012a|[2m3l l. 

7.2 Observations 

7.2.1 Size of the Vela Pulsar’s Radio Emission Region: T = 18 cm 

The fundamental observable of interferometry is visibility, the product of electric fields at a 
pair of antennas ( [Thompson et al.|2001[ l. Because electric fields of all astrophysical sources 
are noiselike, this product must be averaged over some range of time and frequency. For a 
scintillating source, this averaging must be less than the scales of variation of the scintillation 
pattern with time and frequency, to preserve the variation of visibility from scintillation 
( [Gwinn et al.|2000l l. 

For a scintillating point source, in the absence of noise, the distribution of interferomet¬ 
ric visibility is sharply peaked at the origin ( [Gwinn|2001^ . The effect of a small but finite 
emission size is to soften the sharp peak, shift it from the origin, and narrow the distribution. 
As compared with a point-source model, the finite-size distribution without noise peaks at 
larger real part, but has lower probability density at large and small visibility, for the same 
average flux density (or equivalently, the same mean visibility). 

Noise broadens the distribution of visibility. Although noise blurs the distributions and 
their projections, the difference of point-source and finite-size distributions persists, with a 
characteristic W-shaped signature, as Figure[T3]shows . To compare with pulsar observations, 
we must also incorporate the effects of intrinsic variability. Rapid variability modifies the 
noise statistics, while variability over longer times broadens the distribution ( [Gwinn et al.[ 
[201 l[[2012a[l. Consequences of these effects differ from those of emission size. 
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Real Part of Visibility Re[V] Pulse Phase (msec) 


Fig. 13 Observed distribution of visibility 
projected onto the real axis f*, in bins along 
the real axis. Upper: Observed distribution. 
Lower: Residuals to best-fitting model with 
zero size for the Vela pulsar at d = 18 cm. 
Curve shows difference of finite- and zero- 
size models. After|Gwinn et al.]j2Q12b). 


Fig. 14 Best-fitting amplitude (top panel) and source size 
(kMdcr) (lower panel) plotted with pulse gate, for 4 gates in 
6 spectral ranges. The model for the emission region assumes 
a circular Gaussian distribution of emission. After IGwinn eti 

|a:i t2012b) . 


Because finite size narrows the distribution of visibility, and noise broadens it, the dif¬ 
ference of best-fitting models with finite size and zero size has a characteristic W-shaped 
signature. Figure[T3]shows one example of a fit for a range early in the pulse. The character¬ 
istic W-shaped residual is evident, indicating the presence of a finite emission size. A model 
including one additional parameter, for finite size for the pulsar’s emission region, matches 
this residual accurately with significance exceeding 40cr. The inferred size of the emission 
region is 420 km. From fits to gates as a function of pulse phase, we find that the size of the 
pulsar emission region is large at the beginning of the pulse, declines to near zero size near 
the middle of the pulse, and then increases again to nearly 1000 km at the end of the pulse. 
The quoted sizes indicate the full width at half maximum of an equivalent circular Gaussian 
model. 

Theoretical models of pulsar emission typically take their starting point in the geometri¬ 
cal models described above . [Hakoby an & Beskin|p014[ l made theoretical calculations of the 
images of pulsars as a function of pulse phase, using generic expressions for the emission 
altitude and beam shape. They include effects of refraction by the magnetospheric plasma, 
and investigate emission heights up to lOOx the radius of the neutron star. They find a char¬ 
acteristic U-shaped curve of the form seen in Figure [T^ This form results from the greater 
curvature of field lines further from the magnetic pole, and the consequently greater set of 
loci that can emit in a given direction. Interestingly, they find that the size of the emission 
region is much larger than its shift over the course of a pulse. [^en & Melrose| l |2014| l inves¬ 
tigate a similar model, and find that the shift of the emission region over a pulse is indeed 
small. They suggest from geometrical arguments that emission arises at altitudes of more 
than 10% of the light-cylinder radius. [Lyutikov et al.| ( [199^ comes to similar conclusions 
based on emission physics. 
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Fig. 15 Observed and model PDFs of inten¬ 
sity for the subsets of pulses in the top and 
bottom decile by pulse intensity. Theoretical 
residual curves are completely determined 
by a single parameter, the source size. Plot¬ 
ted results correspond to a point source and a 
source that extends over 20 km; the latter is 
clearly inconsistent with the observed statis¬ 
tics. From |Johnson et aT1j2012by 


Fig. 16 Infen'ed emission sizes of individual pulses 
at i = 40 cm. The size is displayed as a function of 
the single-pulse signal-to-noise ratios S ; both linear 
polarizations are plotted. Because S determines the 
standard error for each measurement to excellent ac¬ 
curacy, we omit error bars and instead show the ex¬ 
pected ±3ct errors about = 0. We do not obtain a 
statistically significant detection of emission size for 
any pulse. From|Johnson et ar]j2012by 


7.2.2 Size of the Vela Pulsar’s Emission Region at A = AQ cm from Nyquist-Limited 
Statistics 

The unique nature of pulsar emission allows an elegant solution to determination of the dis¬ 
tribution of intensity of a variable, scintillating source: the formation of spectra that contain 
all single-pulse power. Such spectra require a Fourier transform of a data stream that spans 
the entire pulse, including any scatter-broadening. Without any averaging at all (that is, at the 
Nyquist limit of the data stream), such spectra show the influence of finite source size. |Johm| 
|son & Gwinn| ( [2012a[ l calculated the distribution of Nyquist-sampled spectra, for scintillat¬ 
ing sources with and without effects of size, including the effects of averaging and temporal 
decorrelation. With knowledge of background noise from off-pulse spectra, the intensities of 
individual pulses, and the scintillation timescale, these statistics provide a measure of source 
size. A great strength of this technique is that it can measure size for individual pulses, or 
narrow classes of pulses. 

[Johnson et al.| ( [2012b[ l used the Nyquist-sampled technique to find the size of the Vela 
pulsar at 40-cm wavelength, using baseband recording of the pulsar’s electric field, at the 
Green Bank Telescope. They found that the size was consistent with a pointlike source in 
all cases. The observational upper limit depended upon the set or subset of pulses analyzed. 
Figure [T^ shows a typical example, the distribution of intensity for the brightest 10% of 
pulses, and for the weakest 10%. The distributions are normalized to the mean intensity 
in both cases, so differences arise from the difference in signal-to-noise ratio. The size is 
expressed in terms of the characteristic scales of interstellar scattering by the parameter 
ys = (InMOcrIXf. Here, 6 is the angular broadening by interstellar scatter, and cr is the 
size of a model Gaussian distribution of intensity at the source, both expressed as standard 
deviation. The observing wavelength is A, and M is the ratio of the distance of the observer 
from the scattering screen, to that of the pulsar from the screen. As the figure shows, both 
distributions are clearly inconsistent with a size as large as cTc = 20 km, corresponding to 
a full-width at half-maximum of 47 km of an assumed circular Gaussian emission region. 
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Fig. 17 Left panel: Inferred image of the speckles that scatter pulsar B0834+06. Color is proportional to time 
delay, in a periodic hue map. This image becomes the celestial interferometer aperture, for imaging pulsar 
B0834+06. Center: Pulse profile, with 3 bins indicated. Right: Infened shift of emission region with pulse 
phase. This image from|Pen et al.|j2014) is reproduced couifesy of U.-L. Pen. 


Figure [T^ shows the size of the pulsar as measured in individual pulses, with a range of 
signal-to-noise ratios. From an fit to their full sample of pulses, they obtained a 3cr upper 
limit of cTc < 4 km (FWHM< 9 km). These sizes are comparable to the size of the neutron 
star, and suggest a very concentrated emission region. Theory would suggest that the shift 
of the emission region is still smaller l |Hakobyan & Beskin|2014[|Yuen & Melrose|2014[ l. 

At face value, our results for the size of the emission region of the Vela pulsar at wave¬ 
lengths of 18 and 40 cm appear inconsistent ( |Gwinn et al.||2012b{ [Johnson et al.|2012b^ ). 
How can the size of the emission region change by an order of magnitude, with a change 
of only 2x in observing wavelength? Longer wavelengths are thought to arise at higher al¬ 
titudes, apparently exacerbating the discrepancy ( |Sturrock|1971[ |Ruderman & Sutherland! 
|1975| l. Repeated observations have conhrmed the observational results. 

Refraction of a emergent double-peaked component at 4 = 18 cm may be responsible. 
At 4 = 40 cm the pulse prohle contains a single “core” component, but at 4 = 18 cm 
an additional, double, “cone” component appears ( |Komesaroff et al.|1974[|Kern et al.|2000[ 
[Johnson et al.|2012b^ . As discussed in Section [63| a double-peaked component indicates the 
presence of the 0-mode, and the effects of refraction; whereas a single-peaked component 
indicates the X-mode and no refraction, and consequently a smaller size. Magnetospheric 
refraction might be stronger at the shorter wavelength Arons & Barnard Barnard 


|& Arons] ( |1986[ l; [Lyutikov & Parikh] l |2000[ l; [Hirano & Gwinn (20011; Hakobyan & Beskin 
12014 1 . This matches the observed pattern. 


7.2.5 Femloarcsecond Astrometry of Pulsar B0834+06 

[Pen et af] ( |2014[ l extended the comparison of pulsar scattering patterns at different pulse 
phases. They model the scattering as the interference of a set of points at the screen, the 
“speckles”. Their method isolates the wavefields of each pair of interfering speckles. Inter¬ 
ference of each pair of speckles acts as a 2-slit interferometer to cause the pulsar intensity 
observed at Earth to vary with a specific timescale and bandwidth. A shift of the position 
of the emission region of the pulsar, over the course of a pulse, causes a reflex shift of the 
interference pattern from each pair. Pen et al. inverted the very-long baseline interferometry 
observations of this pulsar by |Brisken et al.| ( |2010[ lto infer the structure of the speckles at 
the scattering screen, as shown in Figure [iTlThey use a holographic technique to partially 
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descatter the data (see also [Walker et al.| l |2008[ l). Their technique corrects for the phase of 
each speckle relative to its neighbors, and so effectively concentrates the power and boosts 
the signal-to-noise ratio. The application of this technique to PSR 0834+06 yields an as¬ 
trometric determination of the phase shift across the pulse profile equivalent to an angular 
resolution of 150picoarcseconds, or 10 km at the distance of the pulsar. This remarkable 
accuracy is comparable to the shift in position of the pulsar due to proper motion, over a sin¬ 
gle pulse. In particular, they found that the velocity of the radio imag e in the picture plane 
is about 1000 km s“’, in good agreement with theoretical prediction (Hakobyan & Beskin 
|M4l l. 


8 Magnetic Axis Alignment 

8.1 Theoretical Predictions for Motion of the Magnetic Axis 

As was shown above, resulting from the MHD theory of the neutron star magnetosphere, 
the magnetospheric torque acting on the surface of a neutron star gives the positive factor 
\Kj_-K^^ in (|12|l-([T3ll corresponding to the alignment evolution of the inclination angle. On 
the other hand, according to this implies that the antisymmetric current fa is to be large 
enough. E.g., for orthogonal rotator the longitudinal current j is to be 10^-10'* times larger 
than the local Goldreich-Julian one yjjf = |fiB|/2;r. Recent simulations of pair production in 
the inner gap ( [Timokhin & Arons|2013[ l suggest that the microphysics of the cascade near 
the polar cap can support the large currents {j » ) required by the global magnetospheric 

structure (it could be accompanied by an efficient heating of the polar cap). Similar results 
are obtained also by global, 3D PIC simulations of pulsar magnetospheres [ Philippov et| 
[ar| ( [2014} . In fact, force-free, MHD, and PIC simulations all find that even though for an 
orthogonal pulsar jgj essentially vanishes due to the midplane symmetry, the magnitude 
of the current flowing out along the open magnetic field lines is very similar to that of the 
aligned pulsar. The results of force-free and MHD simulations tell us that ( [Philippov et al.[ 
[20T4l l, 

^11 — “^aligned COS^, (25) 

Ki^ — ~2A"aligned siu^, (26) 

where ^aligned = Waiigned/f2 = I is the spindown torque of an aligned rotator. There¬ 
fore, 

1,0 ~ -^raligned(l + X)^ (27) 

IrQx ~ -^^aiigned sin^ cos;^- (28) 

Thus, the force-free and MHD simulation results suggest that pulsars tend to become aligned 
with time. This is not surprising in the context of previous discussion: pulsars tend to evolve 
toward the lowest luminosity state, e.g., toward the aligned state (see ( [19[ l. Vacuum pulsars 
become aligned exponentially fast, even before they have a chance to spin down substan¬ 
tially, and generically end up with a period that is a few times their birth period. If most 
pulsars were born as millisecond rotators, this presents a problem in pulsar population syn¬ 
thesis studies, as this would imply that most pulsars would have millisecond periods, yet 
we observe many pulsars with periods of ~second. In contrast, plasma-filled pulsars come 
into alignment much slower, as a power-law in time, x oc so both the spindown and 
alignment proceed at a similar rate ([Philippov et al.|2014[l. 
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On the other hand, if there is some restriction of the value of the longitudinal current 
flowing through the polar cap (no numerical simulation has such a restriction), the situation 
can be different. Such an alternative model in which both symmetric and antisymmetric cur¬ 
rents correspond to the local Goldreich-Julian value was considered by |Beskin et al.| ( [T983| 
|1993| l. They calculated the torque associated with the Ampere force arising from the in¬ 
teraction of the neutron star poloidal field with the surface currents (these currents close 
the longitudinal currents flowing in the region of open magnetosphere). One can prove by 
straightforward but cumbersome calculation that the two approaches are identical. This is 
a crucial assumption because pulsar spindown luminosity is proportional to the magneto- 
spheric current squared. 

For fa » 4 ~ 1 equations (|12|l-(|13|l can be rewritten as 

I,Q ~ cos^x, 

I,Qx ~ /fjj'sin;tcos;t-; 

for orthogonal rotator cos;t' < {QRjcyi^ we have 

40 (31) 

As for cos;t^ > (OJ?/c)’^^ evolutionary equations have an integral 

O sin;^^ = const, (32) 

this model predicts the evolution of the inclination angle toward an orthogonal configura¬ 
tion. Thus, two theoretical models of the neutron star evolution give approximately identical 
predictions for the period derivative P, but opposite ones for the evolution of the inclination 
angle 


(29) 

(30) 


3.2 Observational Constraints on Evolution of Inclination Angle 


Measurement of the rate of change of position angle with pulse phase at the center of the 
pulse dif/ldcf) yields only a measure of the minimum angle between the line of sight and 
the magnetic axis, 4, as inspection of Eqn. \2A\ shows. The angle 4 is sometimes called 
the “impact angle” (see Figure [TOll. Consequently estimates of the inclination angle x are 
indirect, and observational tests of the theories for evolution of^ in the previous section are 
difficult. 

As inspection of Eqn. shows, measurement of the rate of change of position angle 
with pulse phase at the center of the pulse dij/ld^ yields only a measure of the minimum 
angle between the line of sight and the magnetic axis, f, sometimes called the “impact 
angle” (see Figure [T}. Consequently estimates of the inclination angle are indirect, and 
observational tests of the theories for evolution of x in the previous section are difficult. 

In a careful study, |Tauris & Manchester! I |1998[ l compared the inclination angles x and 
rotation period P for nearly 100 pulsars. They found that^ decreases as P increases, over this 
sample of the pulsar population. They made the straightforward assumption that the beam 
from the pulsar is round, as the hollow-cone model discussed in Secti on |2.3| a nd Figure [T] 
suggest]^ They used beam radii as a function of pulse period derived by Gould ( |l994| l from 


■ |Narayan‘& 


Vivekanand 


1983 


decreases as the pulsar ages. 


suggest that pulsar beams are, instead, elongated; and that their elongation 
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Fig. 18 Dependence of the mean inclination angle Fig. 19 Pulsar extinction line in a P- sin^ diagram 
X assi function of the pulsar dynamical age td- Sta- for different magnetic fields. Aitows show the evo- 
tistically this angle decreases with P and, hence, lution tracks of individual pulsars in the model of 
with the age td- After |Tauris & ManchesterH1998), the current losses j29t -j30| jBeskin et al.|1993) . 
using data of|RankinH1993) and|Manchester et al.| 

< 20 ^ . 


comparisons among pulsars with similar periods but different impact angles f, and from pul¬ 
sars with an interpulse (assumed to be nearly orthogonal: x = ^/2) by|Rankin|( |1990| >. From 
the observed pulse width, [Thuris & Manchester| ( [1998| > then inferred the angular separation 
of the line of sight and the rotation axis, and so the inclination angle x- Figure illus¬ 
trates their results, using data from|R ankin| (| 1993] l and |Manchester et al.| (|2005|l.| Weltevrede| 
|& Johnston] ( |2008[ l reached similar conclusions, by comparing the sample of pulsars with 
interpulses with the full population. 


As FigurejT^shows, observations reveal average statistical inclination angles < ;^' > in¬ 
disputably decrease as the period P of pulsars increases and its derivative P decreases. There¬ 
fore, the average inclination angle decreases as the dynamic age increases. Correspondingly, 
pulsars with longer periods exhibit relative ly larger pulse widths Wr = / sin y, where 


VT® is the width of the directivity pattern ( Rankin|l990 


Gould 


1994 


Young et al. 


2010 I. 


These results definitely speak in favor of the alignment mechanism. On the other hand, re¬ 
cently |Lyne et 5^ ( |2013| l on the analysis of the 45 years observations of the Crab pulsar 
concluded that its inclination angle increases with time. Flowever, the effects of stellar non¬ 
sphericity, leading to free precession, can account for this seemingly odd behavior (| Arza-| 
[masskiy et al.|2015b[ l. 

The average inclination angle (x)to for ^ given range of ages Tq can decrease, even if 
the inclination angles of individual pulsars increases with time, in accord with Eqn. ( |32| i. 
This is a consequence of the dependence of the magnetospheric charge density on the incli¬ 
nation angleX- For example, in the picture of [Ruderman & Sutherland| ( |1975| l, radio emis¬ 
sion results from a secondary electron-position cascade, initiated by pair-production from 
curvature-radiation photons. The acceleration of primary electrons within a gap produced 
these curvature photons. The Goldreich-Julian charge density poj ~ QBcosx/( 27 tc) sets the 
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accelerating potential across the gap. As a pulsar ages, pqj and the accelerating potential 
decline with Q. The area of the polar cap also decreases as fewer field lines penetrate the 
light cylinder, and those remaining within the polar cap are less curved. As Eqn. 0 shows, 
the mean free path to pair production increases; the gap becomes wider. When the gap width 
is comparable to the polar-cap radius, the cascade, and radio emission, terminate. 

Because the charge density depends on cos;^^ as well as Bq and Q, death comes to pulsars 
with different inclinations, but the same magnetic field, at different spin periods P. Evalu¬ 
ation of the Ruderman & Sutherland model yields Pmax Indeed, as can 

be seen from Figure [T^ for given values of the pulsar period P and the magnetic field Bq, 
the production of particles is suppressed precisely at inclinations close to njl, where the 
magnetic dipole is nearly orthogonal. Therefore, neutron stars above and to the right of the 
extinction lines in Figure[^do not appear as radio pulsars. 

Because of this dependence of the pulsar extinction line on x, the average inclination 
angles of the observed populations can decrease as the dynamic age increases. A detailed 
analysis, carried out in [Beskin et al.|(|1984|l (see also [Beskin & Nokhrina|(|2004|l;|Beskin| 
|& Elisee^ ( |2005[ )) on the basis of a kinetic equation describing the distribution of pulsars 
confirms this picture quantitatively. 

Clearly, in any case, that the inclination angle is a key hidden parameter: without tak¬ 
ing it into account, it is impossible to construct a consistent theory of the evolution of radio 
pulsars. [Eliseeva et al.j ( |2006 [l include this possibility in their work, suggesting a possible 
direction for further improvements in models for the evolution of neutron stars ([Lipunov et[ 
|al.|1996(|Story et al.|2007[|Popov & Prokhorov|2007[|Gull6n et al.[2014j l. 


9 Summary 

Pulsars provide elegant, although not simple, laboratories for fundamental electromagnetic 
processes at high energies. The basic picture of their structures, involving strong magnetic 
fields and rapid rotation, generation of electron-positron pairs and an energetic wind that 
carries away the rotational kinetic energy of the pulsar, became clear not long after their 
discovery a half-century ago. Recent work has begun to uncover the detailed structures of 
their magnetospheres; the location, size, and optics of their radio emission regions; and evo¬ 
lution of their spins. Some important questions not far from solution include the effect of 
plasma on the magnetosphere and the possible existence of multiple states; conversion of 
Poynting flux to a particle wind (the cr-problem); the location, size, and properties of emis¬ 
sion of different pulse components; and whether rotation and magnetic axes tend to co-align 
or mis-align with time. Further insightful theory, careful statistical studies and targeted ob¬ 
servations will lead to deeper understanding, as pulsars continue their role as the archetypal 
observable neutron stars. 
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